ON 

o\ 



^1- 



o> 

ON 



Clump Giant Distance to the Magellanic Clouds and Anomalous 

Colors in the Galactic Bulge 

Piotr Popowski 

Institute of Geophysics and Planetary Physics, L-413 

Lawrence Livermore National Laboratory, University of California 

P.O.Box 808, Livermore, CA 94551, USA. 



ON I E-mail: popowski@igpp.llnl.gov. 

> 

ABSTRACT 

CO 
CN 



I demonstrate that the two unexpected results in the local Universe: 1) 
anomalous intrinsic (V — I)o colors of the clump giants and RR Lyrae stars 
in the Galactic center, and 2) very short distances to the Magellanic Clouds 



(LMC, SMC) as inferred from clump giants, are connected with each other. The 
(V — I) a anomaly is partially resolved by using the photometry from the phase-II 
of the Optical Gravitational Lensing Experiment (OGLE) rather than phase-I. 
The need for V- or I-magnitude-based change in the bulge (V — I) is one option 
to explain the remaining color discrepancy. Such change may originate in a 
coefficient of selective extinction Ay j 'E(V — I) smaller than typically assumed. 
Application of the (V — I)q correction (independent of its source) doubles the 
slope of the absolute magnitude - metallicity relation for clump giants, so that 
M^RC) = -0.23 + 0.19 [Fe/H]. Consequently, the estimates of the clump 
distances to the LMC and SMC are affected. Udalski's (1998c) distance modulus 
^ ■ of hlmc — 18.18 ±0.06 increases to \ilmc — 18.27±0.07. The distance modulus 

to the SMC increases by 0.12 to Hsmc — 18.77 ± 0.08. I argue that a more 
comprehensive assessment of the metallicity effect on Mj(RC) is needed. 

Subject Headings: distance scale -- dust, extinction -- Galaxy: center - 
Magellanic Clouds — stars: horizontal-branch 



1. Introduction 

Most of the extragalactic distance scale is tied to the LMC, and so the distance to the 
LMC (cIlmc) influences the Hubble constant, Hq. For many years now there has been a 
division between the so called "short" and "long" distance scales to the LMC. Currently, 
the measured values of d^MC span a range of over 25% (see e.g., Feast & Catchpole 1997; 
Stanek, Zaritsky, Sz Harris 1998). Paczyhski Sz Stanek (1998) pointed out that red clump 
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giants should constitute an accurate distance indicator. Udalski et al. (1998a) and Stanek 
et al. (1998) applied the clump method and found a very short distance to the LMC 
{.Hlmc ~ 18.1). In response, Cole (1998) and Girardi et al. (1998) suggested that clump 
giants are not standard candles and that their absolute I magnitudes, Mi(RC), depend on 
the metallicity and age of the population. Udalski (1998b, 1998c) countered this criticism 
by showing that the metallicity dependence is at a low level of about 0.1 mag/dex, and that 
the Mj(RC) is approximately constant for cluster ages between 2 and 10 Gyr. Stanek et 
al. (1999) and Udalski (1999) found a moderate slope of the M^RC) - [Fe/H] relation of 
0.15 mag/dex. The only clump determination, which resulted in a truly long d^MC was a 
study of the field around supernova SN 1987A by Romaniello et al. (1999). However, they 
assumed a bright Mj(RC) from theoretical models and, additionally, the use of the vicinity 
of SN 1987A may not be the most fortunate choice (Udalski 1999). 

The value of Mj(RC) in different stellar systems is a major systematic uncertainty in 
the clump method. It is very hard to prove the standard character of a candle's luminosity. 
However, it should be possible to check whether other stellar characteristics of a candle 
behave in a predictable fashion. Therefore, in §2 I discuss the (V — J)o colors of the clump 
giants and RR Lyrae stars in the Galactic bulge. After making photometric corrections, I 
argue that the remaining color discrepancy between the Baade's Window and local stars 
might have been caused by an overestimated coefficient of selective extinction. Using 
corrected colors, in §3 I derive a new Mj(RC) - [Fe/H] relation for red clump stars and 
show its substantial impact on the distances to the Magellanic Clouds. I summarize the 
results in §4. 



2. Mystery of anomalous colors in the Galactic bulge 

Paczyhski (1998) tried to explain why the clump giants in the Baade's Window 
have (V — J)o colors which are approximately 0.2 magnitudes redder than in the solar 
neighborhood (Paczyhski & Stanek 1998). Paczyhski (1998) suggested super-solar 
metallicities of the Galactic bulge stars as a possible solution. However, there is a 
spectroscopic evidence (see Minniti et al. 1995) that the average metallicity of the bulge is 
[Fe/H] e (—0.3,0.0). Stutz, Popowski & Gould (1999) found a corresponding effect for the 
Baade's Window RR Lyrae stars, which have (V — I)q redder by about 0.17 than their local 
counterparts (Fig. 1). 

A similar size of the color shift in RR Lyrae stars and clump giants suggests a common 
origin of this effect. Does there exists any physical mechanism that could be responsible 
for such behavior? The bulge RR Lyrae stars and clump giants both burn Helium in their 
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cores, but the similarities end here. RR Lyrae stars pulsate, clump giants do not. RR Lyrae 
stars are metal-poor, clump giants are metal-rich. RR Lyrae stars are likely to be a part 
of an axisymmetric stellar halo (e.g., Minniti 1996; Alcock et al. 1998a), whereas clump 
giants form a bar (e.g., Stanek et al. 1994; Ng et al. 1996). For RR Lyrae stars, Stutz 
et al. (1999) suggested that their very red (V — 7)o might have resulted from an unusual 
abundance of a- elements. Why should a clump population which emerged in a different 
formation process share the same property? 

The solutions to the anomalous colors proposed by Paczyhski (1998) and Stutz et 
al. (1999) are not impossible but are rather unlikely. Alternatively, the effect might be 
unrelated to the physics of those stars. The investigated bulge RR Lyrae stars and clump 
giants share two things in common. First, photometry of both types of stars comes from 
the OGLE, phase-I, project. Indeed, Paczyhski et al. (1999) showed that the OGLE-I 
V-magnitudes are 0.021 mag fainter and I-magnitudes 0.035 mag brighter than the better 
calibrated OGLE-II magnitudes. Therefore, the correct (V — I) colors should be 0.056 bluer. 
Additionally, the new (V — I)q from the more homogeneous Baade's Window clump is bluer 
than Paczyhski's & Stanek's (1998) color even when reduced to OGLE-I calibration^. When 
the new OGLE-II photometry reported by Paczyhski et al. (1999) is used, the (V — I)o 
anomaly shrinks and the remaining unexplained shift amounts to ~ 0.11 both for the RR 
Lyrae stars and clump giants. 

Second, Paczyhski & Stanek (1998) and Stutz et al. (1999) use the same extinction 
map (Stanek 1996) and the same coefficient of conversion from visual extinction Ay to 
color excess E(V — I). The absolute values of Ays are likely approximately correct (see 
equation P because the zero point of the extinction map was determined from the (V — K) 
color and Ay/E{V — K) is very close to 1 (Gould, Popowski, & Terndrup 1998; Alcock et 
al. 1998b). However, R VI = A V /E(V — I) is not as secure and has a pronounced effect on 
the obtained color. 

Most of the current studies of the Galactic bulge use Ryj = 2.5. If a true Ryj towards 
Baade's Window equals a instead, then the adjusted Stanek's (1996) V-band extinction, 
will beR 

OL 
^4y,adjusted = TTT \Av ~ ^V,0-point) + ^4v,0-point, (1) 

where most of the extinction, namely v4y -point is excluded from the adjustment because it 



1 Udalski's (1998b) data for the LMC, SMC, and Carina galaxy come from OGLE-II and therefore do not 
require any additional adjustment. 

2 Equation (Q) implicitly assumes that differential (V — I) colors from Stanek (1996) are correct. Whether 
it is the case is an open question. 
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has been determined based on i^-magnitudesQ. The adjustment to the color, which follows 
from equation ([[]) is: 

A(V — J)o = -Z~zA V ^4y,adiusted = —^TZ -^V.O-point • (2) 

2.5 a 2.5 a 

Therefore, for a color shift A(V — I)o , one expects: 

2 .5 Ay : Q_ po i nt 

a ~ A Vfi . point - 2.5 A(V - I) ' U 

Using A(V — J)o ~ —0.11 as required to resolve the color conflict in Baade's Window and 
^v,o-point = 1-37 (Gould et al. 1998; Alcock et al. 1998b), I find a « 2.1 (Fig. 2). This 
i?y/ = 2.1 is certainly low, but not unreasonably so. Szomoru & Guhathakurta (1999) find 
that cirrus clouds in the Galaxy have extinctions consistent with Ay/E(B — V)<2, which 
is more extreme than the change suggested here. If the extinction towards Baade's Window 
is in part provided by the cirrus clouds, then the low Ryi would be expected rather than 
surprising. 

The value and variation of Ryi was thoroughly investigated by Wozniak & Stanek 
(1996). The essence of the Wozniak & Stanek (1996) method to determine differential 
extinction is an assumption that regions of the sky with a lower surface density of stars have 
higher extinction. Wozniak & Stanek (1996) used clump giants to convert a certain density 
of stars to an amount of visual extinction. To make a calibration procedure completely 
unbiased would require, among other things, that clump giants were selected without any 
assumption about Ryi', that absolute V^-magnitudes of clump giants, My(RC), do not 
depend on their color [here (V — I)o]\ and that reddened and unreddened clump giants 
be drawn from the same parent population. None of these is true. A color-magnitude 
diagram (CMD) for dense Galactic fields does not allow one to unambiguously distinguish 
clump giants from other stars. Different parts of an intrinsically clean CMD overlap 
due to differential reddening and a range of stellar distances. Therefore, the selection of 
clump giants must involve some assumptions about Ryi- Wozniak & Stanek (1996) adopt 
Rvi = 2.6. This procedure tends to bias the derived relation toward this predefined slope. 
Wozniak & Stanek (1996) were fully aware of this effect, and they performed a number of 



3 To the zero-th order, a shape of the extinction curve is governed by one parameter, so that the change 
in R V i will affect A v / E(V - K). In effect, Gould et al. (1998) and Alcock et al. (1998) might have 
overestimated (underestimated) ^4y : o-point by a few percent if a < 2.5 (a > 2.5). This change would 
propagate to equations (|l|) - (g) in several correlated ways. Thus, strictly speaking a value of a is additionally 
a function of Ay / E(V — K), which changes with a. Therefore, solving this problem exactly requires iteration. 
I neglect this complication in further considerations. 
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simulations, which are summarized in their Figure 4. In brief, in the range 2.1 < Ryj < 3.1, 
the bias scales as SRyi ~ 0.4 (2.6 — Ryi) and so may become very substantial for a low 
or high Ryi- In particular, if the true Ryj = 2.1, Wozniak and Stanek (1996) would find 
Ryi = 2.3. Therefore, this effect alone could account for half of the difference between the 
required and measured Ryi- 

The intrinsic characteristics of the bulge clump stars are unknown, but I will assume 
they resemble the clump measured by Hipparcos (European Space Agency 1997). The fit to 
the local clump giants selected by Paczyhski & Stanek (1998) gives My(RC) oc 0.4(V — 1) . 
Therefore, the structure of the local clump itself acts similarly to extinction with Ryi = 0.4. 
In an ideal case, when the CMD locations of the entire clump populations in different 
fields are compared, the M V (RC) — (V — I) dependence should not matter. However, 
when combined with the actual extinction and additionally influenced by the completeness 
function of a survey, this effect may additionally bias the value of Ryj. 

Because the smaller selective extinction coefficient is not excluded by the current 
studies, one can assume Ryj = 2.1 to match the (V — I)q colors of the bulge with the ones 
in the solar neighborhood. The color is a weak function of [Fe/H], so this procedure is 
justified because the [Fe/H] of the bulge and solar neighborhood are similar. This change 
in Ryi will decrease the I-mag extinction, Ai, by 0.11 mag. Therefore, the clump distance 
to the Galactic center would increase by the same amount. 



3. Recalibration of the clump 

What is the bearing of the bulge results on the distance to the LMC? Let A indicate the 
difference between the mean dereddened I-magnitude of clump giants and the derredened 
V-magnitude of RR Lyrae stars at the metallicity of RR Lyrae stars in the Galactic bulge. 
When monitored in several stellar systems with different clump metallicities, the variable 
A, introduced by Udalski (1998b), allows one to calibrate the Mj(RC) - [Fe/H] relation 
with respect to the baseline provided by RR Lyrae stars. The better photometry from 
Paczynski et al. (1999) and a possible modification of Ryi influence the value of A at 
the Galactic center (Abw)- It is important to note that one will face the same type of 
adjustment to Abw whenever the anomalous colors in the Baade's Window are resolved at 
the expense of the modification of V- or I-magnitudes. That is, the modification of Ryi is 
not a necessary condition! It is simply one of the options. As a result of the change in Abw, 
the Mi(RC) - [Fe/H] relation for clump giants changes. Moreover, \ilmc and hsmc will 
change as well because the Mj(RC) - [Fe/H] relation is used to obtain the clump distances 
to the Magellanic Clouds. 
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Here, I will modify Udalski's (1998b) A versus [Fe/H] plot and derive a new Mj(RC) 
- [Fe/H] relation consistent with the new data and considerations from §2. I construct the 
Udalski (1998b) plot using his original points modified in the following way: 

- To match the change in (V — I)q, I modify Abw by 0.17 mags (a combined change from 
photometry and some other, yet unrecognized, source, e.g., selective extinction coefficient). 

- I modify the [Fe/H] of the Baade's Window clump giants, so that [Fe/H] = 0.0 (see e.g., 
Minniti et al. 1995 for a review on the bulge metallicity). 

The possible improvement to the above procedure would be a construction of Udalski's 
(1998b) diagram based on clump giants in the LMC and SMC clusters, which would reduce 
the uncertainties associated with the reddening to the field stars. This more complex 
treatment is beyond the scope of this paper. 

I make a linear fit to the A - [Fe/H] relation. I assume that a total 
error in dependent variable A for the z-th point, a to tai,ii can be expressed as 

/ \ 2 

.2 _ ^.2 i / dA 



Vtotal,i = <7A,i + ^dlFVHl [Fe/HI^/^'V ' ^^ ^ ^ <T ^>f S i> i ^ ^ illdividual 

point errors in A and metallicity, respectively. As a result, I obtain: 

M T (RC) = (-0.36 ± 0.03) + (0.19 ± 0.05) ([Fe/H] + 0.66) (4) 

Equation (J|) is expressed in the form with uncorrelated errors and normalized to the local 
Hipparcos (European Space Agency 1997) result of Mi(RC) = -0.23 ± 0.02 at [Fe/H] 
= 0.0 reported by Stanek & Garnavich (1998). Equation (f|) is a good fit to the data with 
a x 2 /d-o.f = 3.17/2 (Fig. 3). The slope from equation (f|) is 0.10 mag/dex steeper than 
the one given by Udalski (1998b). However, it agrees well with the slope of 0.15 ± 0.05 
based on the spectroscopic data for the local clump (Udalski 1999). Also, the slope from 
equation (11) is in a good agreement with the theoretical models (e.g., Cole 1998) and leads 
to M I)LMC {RC) = -0.35 in the LMC field and M ItLMC {RC) = -0.39 for the clusters (close 
to the values suggested by Girardi et al. 1998). 



4. Summary 

I demonstrated that the correction of the bulge (V — I)o anomaly has a pronounced 
effect on the slope of the absolute magnitude - metallicity relation for clump 
giants. Introducing color correction to the original Udalski (1998b) diagram, I find 
M^RC) = -0.23 + 0.19 [Fe/H]. Consequently, the Udalski's (1998c) distance modulus of 
Hlmc — 18.18 ± 0.06 is increased to [Ilmc — 18.27 ± 0.07. The distance modulus to the 
SMC increases from /i SMC = 18.65 ± 0.08 to fi S MC = 18.77 ± 0.08. 
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Even though my approach in this paper is only qualitative, there are two important 
characteristics of this study: 

1) The calibration of Mj(RC) - [Fe/H] relation, has been based on the homogeneous set of 
the OGLE-II photometry. Therefore, no corrections due to the use of different telescopes, 
instruments and reduction procedures are required. Unfortunately, this makes the above 
calibration vulnerable to unrecognized systematic problems of the OGLE photometry. 

2) The Mi(RC) value has been derived based on observational data and not simply picked 
from a family of possible theoretical models of stellar evolution. 

Romaniello et al. (1999) provide an independent source of clump photometry in the LMC, 
but due to the importance of photometric homogeneity I am not able to use their data in 
a way consistent with the rest of my analysis. With reference to point 2), it is crucial to 
note that observationally calibrated Mj(RC) is not subject to the modeling uncertainties 
which affect the Romaniello et al. (1999) distance to the LMC. However, my calibration is 
only as good as the assumptions and data that enter the analysis. Reddening corrections 
to the original Udalski's (1998b) diagram, which is partly based on the field stars in the 
LMC and SMC, may be needed. Therefore, a more comprehensive study of the metallicity 
effect on Mi(RC) is necessary. Udalski's (1999) determination based on the local clump is 
an important step toward establishing a reliable Mj(RC) - [Fe/H] relation. 

Andy Becker deserves my special thanks for many stimulating discussions about the 
extinction issues in the Galactic bulge. I am deeply grateful to Andy Gould for his very 
careful reading of the original version of this paper and a number of insightful remarks. 
I also would like to thank Kem Cook for his valuable comments and discussions. Work 
performed at the LLNL is supported by the DOE under contract W7405-ENG-48. 



REFERENCES 

Alcock, C. et al. 1998a, ApJ, 492, 190 

Alcock, C. et al. 1998b, ApJ, 494, 396 

Cole, A.A. 1998, ApJ, 500, L137 

European Space Agency. 1997, The Hipparcos Catalogue (ESA SP-1200) (Paris: ESA) 

Feast, M.W., & Catchpole, R.M. 1997, MNRAS, 286, LI 

Girardi, L., Groenewegen, M.A.T., Weiss, A., & Salaris, M. 1998, MNRAS, 301, 149 

Gould, A., Popowski, P., & Terndrup, D.M. 1998, ApJ, 492, 778 



Minniti, D. 1996, ApJ, 459, 175 

Minniti, D., Olszewski, E.W., Liebert, J., White, S.D.M., Hill, J.M., & Irwin, M.J. 1995, 
MNRAS, 277, 1293 

Ng, Y.K., Bertelli, G., Chiosi, C, & Bressan, A. 1996, A&A, 310, 771 

Paczyhski, B. 1998, Acta Astron., 48, 405 

Paczyhski, B., & Stanek, K.Z. 1998, ApJ, 494, L219 



Paczyfiski, B., Udalski, A., Szymafiski, M., Kubiak, M., Pietrzyhski, G., Soszyrlski, I. 
Wozniak, P., & Zebrurl, K., 1999, Acta Astron., 49, 319 

Romaniello, M., Salaris, M., Cassisi, S., & Panagia, N. 1999, accepted to ApJ flastro- 



pri/9910082j) 



Stanek, K.Z. 1996, ApJ, 460, L37. 

Stanek, K.Z., & Garnavich, P.M. 1998, ApJ, 500, L141 

Stanek, K.Z., Kaluzny, J., Wysocka, A., & Thompson, I. 1999, submitted to AJ 
( |astr"o-ph/990804l|) 

Stanek, K.Z., Mateo, M., Udalski, A., Szymahski, M., Kaluzny, J., & Kubiak, M. 1994, 
ApJ, 429, L7 

Stanek, K.Z., Zaritsky, D., & Harris, J. 1998, ApJ, 500, L141 

Stutz, A., Popowski, P., & Gould, A. 1999, ApJ, 521, 206 

Szomoru, A., & Guhathakurta, P. 1999, A J, 117, 2226 

Udalski, A. 1998b, Acta Astron., 48, 113 

Udalski, A. 1998c, Acta Astron., 48, 383 

Udalski, A. 1999, submitted to ApJ Letters (|astro-ph/9910167|) 



Udalski, A., Szymahski, M., Kubiak, M., Pietrzyhski, G., Wozniak, P., & Zebruh, K. 1998a, 



Acta Astron., 48, 1 



Walker, A.R. 1992, ApJ, 390, L81 
Wozniak, P., & Stanek, K.Z. 1996, 464, 233 



This preprint was prepared with the AAS IATgX macros v4.0. 



-9- 



Fig. 1. — Left part presents (V — I)q colors of the Baade's Window (open circles) and local 
(full squares) RR Lyrae stars as a function of (V — K) . The (V — I)q shift between the 
two groups is of order of 0.17 mag. The standard extinction line (dashed) is parallel to the 
stellar lines and so will not cure the anomaly. The right panel shows a similar situation for 
clump giants. The vertical (V — I)q color axis should be treated as an abscissa. The two 
Gaussians are Paczyhski's & Stanek's (1998) fits to the number densities of clump giants as 
a function of color: in the solar neighborhood (thick solid line) and in the Baade's Window 
(thin dotted line). The (V — I)o difference between the maxima of the two clump groups is 
0.21. 

Fig. 2. — The expected adjustment to the (V — I)o as a function of the selective extinction 
coefficient Ryi- The full dot marks the standard extinction. The dotted lines indicate the 
new point which would explain the entire color anomaly seen in the Baade's Window. 

Fig. 3. — Udalski's (1998b) diagram showing A versus [Fe/H]. In order of increasing [Fe/H], 
the data points correspond to the Carina, SMC, LMC, and Baade's Window. The V- 
magnitude adjustment to RR Lyrae stars is limited to the OGLE-II-based correction in 
Baade's Window. The possible changes to A that move all the points on the diagram by the 
same amount can be ignored because they do not influence the calibration. 
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